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Interstellar Extinction Law toward the Galactic Center IH : 
J, H, Ks bands in the 2MASS and the MKO systems, 
and 3.6, 4.5, 5.8, 8.0 /im in the Spitzer/l'R.AQ system 

Shogo Nishiyama^, Motohide Tamura^, Hirofumi Hatano^, Daisuke Kato^, Toshihiko 
Tanabe^, Koji Sugitani^, and Tetsuya Nagata^ 

ABSTRACT 

We have determined interstellar extinction law toward the Galactic cen- 
ter (GC) at the wavelength from 1.2 to 8.0 yum, using point sources de- 
tected in the IRSF/SIRIUS near-infrared survey and those in the 2MASS and 
Spitzer/YRkC/GUMP^E H catalogs. The central region | Z |< 3?0 and | 6 |< 1?0 
has been surveyed in the J, H and Ks bands with the IRSF telescope and the 
SIRIUS camera whose filters are similar to the Mauna Kea Observatories (MKO) 
near-infrared photometric system. Combined with the GLIMPSE II point source 
catalog, we made Ks versus Ks — A color-magnitude diagrams where A = 3.6, 
4.5, 5.8, and 8.0 yum. The Ks magnitudes of bulge red clump stars and the 
Ks — A colors of red giant branches are used as a tracer of the reddening vector 
in the color-magnitude diagrams. From these magnitudes and colors, we have 
obtained the ratios of total to selective extinction Axg/ E^s-x for the four IRAC 
bands. Combined with Ax/A^s for the J and H bands derived by Nishiyama 
et al., we obtain Aj : Ah : Aks ■ ^[3.6] : ^[4.5] : ^[5.8] : ^[s.o] = 3.02 : 1.73 : 1 : 
0.50 : 0.39 : 0.36 : 0.43 for the line of sight toward the GC. This confirms the 
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flattening of the extinction curve at A > S/im from a simple extrapolation of the 
power-law extinction at shorter wavelengths, in accordance with recent studies. 
The extinction law in the 2MASS JHKs bands has also been calculated, and a 
good agreement with that in the MKO system is found. Thus, it is established 
that the extinction in the wavelength range of J, H, and Ks is well fitted by 
a power law of steep decrease Ax oc A~^ ° toward the GC. In nearby molecular 
clouds and diffuse interstellar medium, the lack of rehable measurements of the 
total to selective extinction ratios hampers unambiguous determination of the 
extinction law; however, observational results toward these lines of sight cannot 
be reconciled with a single extinction law. 

Subject headings: infrared: ISM — dust, extinction — stars: horizontal-branch 
— Galaxy: center 

1. Introduction 

The absolute value of interstellar extinction to an individual star is difficult to determine, 
and so is its wavelength dependence. Toward the Galactic center (GC), however, we can 
directly derive the wavelength dependence of interstellar extinction, assuming only that the 
center of stellar distribution in the lines of sight is at the same distance from us and that the 
foreground extinction is patchy; such a principle has been employed for star clusters and is 
referred to as the cluster method or variable-extinction method (Krelowski & Papaj 1993). 
By plotting the apparent magnitude versus the color excess of a group of stars, one obtains a 
straight line with the slope equal to the total to selective extinction ratio, e.g., Axg/ Eh-Ks- 
This variable-extinction method was applied for the red clump (RC) stars around the GC 
by Wozniak & Stanek (1996) at V and /, and by Nishiyama et al. (2006, 2008) as the "RC 
method" at V , J, H, and Ks- 

Nishiyama et al. (2006) measured the peak position of the RC stars in the color- 
magnitude diagram (CMD) of each small (~ 4' x 4') field in their survey area of | / |< 2° 
and I b |< 1°, and the reddening and the extinction were derived from the relative shifts of 
the peak color and the peak magnitude, respectively. Here the most important point is that 
the magnitude of extinction accompanied by a certain amount of reddening can be derived 
exactly, as long as no difference in 1) the mean distance, 2) the mean absolute magnitude, 
and 3) the mean color of the RC population exists among the small fields. Since the survey 
is deep enough and the survey area of a few degrees (projected area of a few hundred pc at 
the GC) does not seem to have different RC population correlated with the amount of fore- 
ground extinction, this assumption is reasonable. The ratios of total to selective extinction. 
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AksI Eh-Ks-i -^KsI Ej-Ks-i ^iid Ah/Ej-h were in turn used to derive the wavelength depen- 
dence of extinction Aj : Afj : Ak^. When approximated as a power law of the wavelength 
A, the interstellar extinction Ax toward the GC decreases as A~^'^^ in the wavelength range 
of J, H, and Ks, being steeper to the longer wavelength than the Rieke & Lebofsky (1985) 
result A~^"^. 

Direct measurements of extinction are extremely difficult because of the uncertainty of 
the distance of light sources. Therefore, the ratios of reddening (ratios of color excesses; 
e.g., Ex-v/Eb-v) ^^'^ generally determined first, and then these color excess ratios are 
transformed into the absolute extinction Ax with some assumptions. For this purpose, 
the total to selective extinction ratio Av/Eb-v = Rv is often derived somehow, either by 
extrapolation of the color excess diagram to the longer wavelength or by a certain assumption, 
and then employed with the equation 

^^-^ + Rv. (1) 



Eb-v Eb-v 

So long as one compares the color excesses, which are obtained directly and accurately, 
meaningful comparison is possible for different determinations, and many attempts have 
been made in this color-color method. However, caution must be exercised in the use of 
absolute extinction ^4^. 

Rickc & Lebofsky (1985) determined the extinction law toward the GC from 1 to 13 
lj,m by the color excess observations of five supergiants near the GC and o Sco. They first 
set Ev-k/Eb-v — 2.744 for all the stars and assumed Ay /Eb-v — 3.09 ± 0.3 for the 
extinction toward the GC. The lower limit of total to selective extinction ratio 3.09 was 
determined from the decrease of extinction in the range of L, M, 8 /im, but the observation 
at these wavelengths was relatively uncertain. Longward of 3 /xm, Lutz ct al. (1996) and 
Lutz (1999) observed hydrogen recombination line emission toward the GC with the Infrared 
Space Observatory SWS instrument and found that the extinction curve is much fiatter than 
the Rieke & Lebofsky (1985) results. Lutz et al. (1996) compared the observed fiuxes and 
expected fluxes predicted for the HE region in an aperture of 14" x 20" centered on Sgr A*. 
They confirmed the applicability of the case B conditions from the observations of different 
upper-level lines. It should be noted, however, that Lutz et al. (1996) had to assume a K 
band extinction of 3.47 mag because the ratio of reddening, not the absolute extinction, in 
the hydrogen line strengths were derived also in their determination. 

The wavelength dependence of interstellar extinction toward the GC, therefore, can be 
summarized as follows. The extinction in the near-infrared wavebands J, H, Ks is fitted well 
by a steep power law X"^, but decreases only slightly beyond A = 3 /im, and then shows a 
large maximum due to the silicate absorption at A = 9.7 /im. The steep power law in the 



-4- 



wavelength range 1 — 2.5 //m and flat curve in 3 — 4 //m is consistent with the polarimetry 
results (Nagata et al. 1994), which can be regarded, in a sense, as an absolute determination 
of the difference between the extinctions of two orthogonal axes. 

In directions of other than the GC, flat extinction dependence in the wavelength range 
of 3 — 8 fjLm has also been derived. Indebetouw et al. (2005) measured the mean color 

excess ratios from the color distributions of stars observed with the Infrared Array Camera 
(IRAC; Fazio et al. 2004) on board Spitzer Space Telescope (SST), along two lines of sight 
in the Galactic plane. Investigations toward five nearby star-forming regions (Flaherty et al. 
2007) and the star-forming dense core Barnard 59 (Roman-Zuhiga et al. 2007) also suggest 
relatively flat extinction curves from 3 to 8 /xm. All these studies obtained color ratios, either 
E\-Ks/ Ej-Ks or Ex-Ks/^H-Ksj first. Then these color excess ratios are transformed into 
the absolute extinction ratio (e.g., Ax/Axg) on the assumption of one ratio (e.g., Ah/Aks) 
using the equation 



Assuming their distribution in the Galactic plane, Indebetouw et al. (2005) fitted the locus 

of RC stars in a color magnitude diagram and derived Ah/A^s = 1-55 ± 0.1. Then using 
the equation (2), Indebetouw et al. (2005), Roman-Zufiiga et al. (2007), and Flaherty et al. 
(2007) derived the absolute extinction ratios. As pointed out by Indebetouw et al. (2005), 
such extinction ratios of Ah/Aks on the assumption of the RC star location in the Galactic 
plane might be more accurate than the extrapolation of the Ex-Ks/ Ej-Ks curve toward the 
longer wavelength to get the absolute extinction values as in equation (1), but one should 
keep in mind that the assumption of the one ratio can lead to large errors in Ax/A^s- Color 
excess ratios with the 5'5'r/IRAC bands are different between the star forming regions and 
the diffuse interstellar medium, and a difference in the extinction law between them was 
suggested (Flaherty et al. 2007). 

In this paper, we assume that the center of distribution in the lines of sight is at the 
same distance from us for the RC giants and the giants in the upper red giant branch (RGB), 
and determined the total to selective extinction ratios A^^/ Ek^-x for the IRAC wavebands 
while we use the IRSF/SIRIUS (similar to the the Mauna Kea Observatories (MKO) near- 
infrared photometric system; Tokunaga et al. 2002) survey results for the Kg band. Also, 
using the RC and upper RGB stars toward the GC, the total to selective extinction ratios 

^Ks/Ej^MAss-Ks, ^Ks/Eh^mass-Ks, ^ud Aks / Ek^^^ss-Ks 1^ thc Two Micron All Sky Sur- 
vey (2MASS) bands were derived, and the agreement to the steep extinction law determined 
by Nishiyama et al. (2006) is examined in the 2MASS bands. We compare these results with 
the previous determination of interstellar reddening in different lines of sight, and discuss 
their difference. 




(2) 
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This variation of the "RC method" is unique and different from any previous determi- 
nations of mid- infrared extinction because it is free from the transformation equations (1) 
and (2). It is on the assumptions of the same mean distance and magnitude of the RC star 
population among the small fields, of the agreement of the central positions of the spatial 
distributions of the RC and RGB stars, and that the upper RGB stars have colors varying 
only with luminosity. The distance to the RC stars seem to be rather constant (Nishiyama 
et al. 2005) in this range of Galactic longitudes, and it seems reasonable to assume that 
the RC and RGB stars have similar spatial distributions. The upper RGB has been used 
for determining the interstellar reddening and extinction by many investigators. Progel et 
al. (1999) used RGB stars with the unreddened K magnitude range of 8.0 to 12.5. They 
compared the RGB colors of the inner bulge (|/| < 4°, |6| < 3°) with that of the Baade's 
Window (Ticdc et al. 1995), and derived the interstellar extinction. They also found that 
the amplitude of metallicity variations in the inner bulge does not cause large RGB slope 
changes in the J — K versus K CMD. Schultheis et al. (1999) used calculated isochrones as 
reference in drawing the extinction map of the inner bulge, and Dutra et al. (2003) again 
used the Baade's Window data as reference in their determining the extinction within 10° 
of the GC. Note that although we have no reliable reference of reddening-free RGB in the 
SST /IRAC bands, only the relative shifts of magnitude and color are needed in the current 
work. 



2. Observational Data 

2.1. IRSF/SIRIUS 

The central region of our Galaxy, | / |< 3?0 and | 6 |< 1?0 (Fig. 1), was observed from 
2002 to 2004 using the NIR camera SIRIUS (Simultaneous Infrared Imager for Unbiased 
Survey; Nagashima et al. 1999; Nagayama et al. 2003) on the 1.4 m telescope IRSF (Infrared 
Survey Facility). The SIRIUS camera can provide J (1.25 /im), H (1.63 /im), and Ks (2.14 
lim) images simultaneously, with a field of view of 7^7 x 7^7 and a pixel scale of 0'.'45. 

Data reduction was carried out with the IRAF (Imaging Reduction and Analysis Facil- 
ity)^ software package. Images were pre-reduced following the standard procedures of near- 
infrared arrays (dark frame subtraction, flat fielding, and sky subtraction). Photometry, 
including point spread function (PSF) fitting, was carried out with the DAOPHOT package 



^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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(Stetson 1987). We used the DAOFIND task to identify point sources, and the sources were 
then input for PSF-fitting photometry to the ALLSTAR task. About 20 sources were used 
to construct the PSF for each image. 

We performed a photometric cahbration with 2MASS point source catalog (Skrutskie 
et al. 2006). When fitted with the Gauss function, the histograms of the difference of the 
2MASS and SIRIUS magnitudes (m2MAss ~ "^sirius) in the three bands for more than 10^ 

stars have standard deviations of ^ 0.03, suggesting that the accuracy of the zero-point 
cahbration for each star is about 0.03 mag. The averages of the lOcr hmiting magnitudes 
were J=17.1, H = 16.6 and Ks = 15.6. 

Indication of the internal reliability of our photometry is obtained from overlapped 
regions between adjacent fields. In our observations, images were taken under different sky 
conditions and at different nights, even in different years. The variation in photometry was 
found due to the different PSF models and zero-point correction used for the analysis of 
each field. We have thus calculated magnitude differences of the same stars in the adjacent 
fields. When we use stars whose photometric errors calculated with IRAF are < 0.01 mag, 
the mean and standard deviation of the magnitude difference are less than 0.01 mag and 
0.03—0.04 mag, respectively, in the three bands. 

2.2. GLIMPSE n 

The Galactic bulge region was observed as one of the SST Legacy programs. Galactic 
Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE) II. The region was imaged 
with three 1.2 s exposures at each location using the IRAC, which is a four channel camera 
operating simultaneously in wave bands, [3.6], [4.5], [5.8], and [8.0], centered on 3.6, 4.4, 5.7, 
and 7.9 /im, respectively. 

Two catalogs are provided by the GLIMPSE n project, one is a highly reliable Point 
Source Catalog (GLMIIC), and the other is a more complete Point Source Archive (GLMI 
lA). In this study, we use the GLMIIC. The criteria to be included in the GLMIIC are 
described in Meade ct al. (2008): e.g., detected at least twice in one band, at least once in 
an adjacent band. The 5 a limiting magnitudes of the point sources are approximately 14, 
12, 10.5, and 9.0 mag in the [3.6], [4.5], [5.8], and [8.0] bands, respectively. Since the central 
region ( | / |< 1° and | b |< 0?75 ) was observed in another program, the list of sources in 
this region is not included in the GLMIIC. 
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3. Reduction and Analysis 

The stars found with IRSF/SIRIUS and in the GLMIIC have been cross-identified using 
a simple positional correlation. The identification was performed with a search radius of 
0'/6, and we found ~ 5.3 x 10^ matches with an rms error less than 0'.'2 in R.A. and Dec, in 
the difference between SIRIUS and the GLMIIC coordinates. The distribution of matched 
sources is shown in Fig. 1. 




Fig. 1. — A source density map for stars matched between the IRSF/SIRIUS catalog and the 
GLMnC. The bins are 0?032 x 0?032, and star count density is given in unit of [lOOOdeg"^)]. 



Fig. 2 shows one of the star counts of 0?2 x 0?2 fields for IRSF/SIRIUS, 2MASS, and 
SST/IKAC Clear peaks of RC stars are seen in the star counts of IRSF/SIRIUS (indicated 
by arrows), but not seen in those of 2MASS and SST/IKAC. From the limiting magnitudes 
obtained by GLIMPSE I (Benjamin et al. 2003), we expected that RC peaks would be 
detected at least in the [3.6] and [4.5] bands, but source confusion might affect the detection 
limit of these bands. 

To make CMDs involving the SST/IKAC bands when we use the "RC method" following 
Nishiyama et al. (2006), we need to find stars that are detected in the SST/IKAC bands. 
Here, we use the upper RGB stars to derive the reddening. We can use the IRSF/SIRIUS 
Ks magnitudes of the RC stars for the derivation of extinction. Assuming that the RC stars 
and the upper RGB stars are similarly distributed in space with their centers at the GC, 
we will be able to determine in the CMD the extinction (from RC magnitude shift) and the 
reddening (from RGB color shift) of a line of sight toward the GC. By plotting the apparent 
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magnitude of the RC stars and the color excess of the RGB of each small field, we obtain 
a straight line with the slope equal to the total to selective extinction ratio A^g/ E^g-x for 
the SST/IRAC bands. 

As a first step in our analysis, we divided the survey area (Fig. 1) into fields of 0?2 x 0?2. 
We then made a Xg-band luminosity function (LF) for each field. A sample of the LFs is 
shown in Fig. 3, top left panel. A clear peak of RC stars can be found at Ks ~ 14 mag. 
We determined the Ks peak magnitude {Ksrc) by fitting with the Gaussian function (thick 
line on the LF in Fig. 3). 

Second, to determine the color of the RGB, we made Ks versus Ks — A CMDs for each 
field and IRAC band. CMDs of one of the fields are shown in Fig. 4 (see also Fig. 3, top 
right panel for A = 4.5 fim). Third, we determined the RGB slope for each field and each 
IRAC band on the assumption that the RGB can be fitted by a linear function. To do this, 
we divided Ks magnitude of stars on the RGB into bins of equal (0.5 mag) size, which is 
represented in the top right panel. Fig. 3 by dashed rectangles. Then we calculated the 
mean of Ks — A color of stars in each bin, by fitting the color histogram with the Gaussian 
function (the bottom left panel in Fig. 3) . The bins near the limiting magnitudes and those 
without enough number of stars were excluded. The mean Ks — [4.5] colors for the bins 
are shown by open circles in the CMD (the top right panel). The RGB slope for each field 
was obtained from the least-squares fit of the mean Ks — A colors with two free parameters: 
slope and intercept (thick line in the top right panel). The magnitude ranges seem to be 
wide enough to determine the fairly linear part of the RGB, and narrow enough to exclude 
both the brighter asymptotic giant branch stars and the fainter RC stars, in these CMDs^. 
We did not find a correlation between the RGB slopes and extinction for each field. 

Fourth, we calculated the mean of the RGB slopes in all the IRAC bands. The histogram 
of the (slopes)"^ in the Ks versus Ks — [4.5] CMDs is shown in the bottom right panel. The 
means of the RGB slopes were then obtained by fitting the histograms with the Gaussian 
function (the bottom right panel in Fig. 3), as 33.6, 17.9, —24.6, and —19.1, for the [3.6], 
[4.5], [5.8], and [8.0] bands, respectively. 

Fifth, the color of the RGB in each field was determined assuming that the RGBs in all 



^The fitted magnitude range for the [3.6] and [4.5] bands differs sliglitly from field to field, but is ap- 
proximately from {Ksrc ~ 2-7) to {Ksrc ~ 0-8), where Ksrc is the peak i^^-magnitude of RC stars for 
each field, obtained in the first step. Those for the [5.8] and [8.0] bands are between {Ksrc ~ 4-2) and 
{Ksrc ^ 2.2), and between {Ksrc ^ 4.2) and {Ksrc ^ 2.6). As shown by the LFs (Fig. 2), the hmiting 
magnitudes in the [5.8] and [8.0] bands are ~ 2 mag brighter than others, and thus the fitting range is also 
brighter by ~ 2 mag. 
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Magnitude 

Fig. 2.— Star counts in the J,H and Kg bands from IRSF/SIRIUS {top) and 2MASS 
{middle), and in the [3.6], [4.5], [5.8], and [8.0] bands from 5'5'r/IRAC for one of the fields 
h = +0?9, — 0?9) towards the GC. Peaks of RC stars can be found only in the top panel, 
and are indicated by arrows. 
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Ks - [4.5] 



Fig. 3. — Ks-haad luminosity function (top left panel) and Ks versus Ks — [4.5] CMD (top 
right panel) for one of the fields (centered at I = — 0?3, 6 = — 0?7). The mean Ks — [4.5] 
colors of stars in the dashed rectangles arc shown by open circles. By fitting the mean colors, 
we can determine the RGB slope for this field. Bottom left panel : Kg — [4.5] histograms of 
stars at 13.0 < Ks < 13.5, which is shown by the bottom rectangle in the top right panel. 
The histogram is fitted with the Gaussian function to obtain the mean of color. Bottom 
right panel : Histogram of (RGB slope)^^ in the Ks versus Ks — [4.5] CMDs of the fields. 
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Fig. 4.— Kg versus Kg-X CMDs of one field (centered at I = -0?3, b = -0?7) for A = [3.6] 
(top left), [4.5] (top right), [5.8] (bottom left), and [8.0] (bottom right). 
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the fields have the same slope as that obtained in the third and fourth steps (the mean of the 
RGB slope). The mean Ks — A color for each field was determined by a one-parameter (i.e., 
intercept) least-squares fit with the result expressed as the color the RGB stars would have 
at the magnitude of the RC stars. Because of the assumption of constant slope, what was 
actually measured is the color at the centroid of the RGB color-magnitude distribution. The 
coordinate (i^5— A, Ksrc) is the indication of reddening and extinction of each field although 
RC stars are undetected in the GLIMPSE n catalog (and also in the 2MASS catalog). 

4. Results 

Fig. 5 (a)-(d) shows the Ks magnitude of the RC peak Ksrc ^"^^ the [Ks — A) color 
of the RGB of each field in Ks versus Ks — A CMDs. Error bars on y-axcs show uncer- 
tainty in RC peak determination, and those on x-axes come from uncertainty of intercepts 
when the least-squares fit is adopted for the mean Kg — A colors. These errors seem to be 
underestimated because the dispersion of the data points is large compared with the error 
bars. Hence we estimated the slopes Axg/Exs-x and their errors by applying a least-squares 
linear fit with minimization under the assumption that the errors are all equal. We show 
the linear fits applied to the data points in Fig. 5, and the resultant slopes and their errors 
are listed in Table 1. 

The same procedure described in the previous section was apphed to 2MASS point 
sources in the J, H and Ks bands, which are included in the GLMIIC. Fig. 6 show the Kg 
magnitude of the RC peak and the relative A2MASS ~ Ks color of the RGB in Ks versus 
A2MASS ~ CMDs. The dependence of -ft'52MASS ~ -^5 '^^ ^s is very small; hence the same 
plot in a CMD was made for Ks and -ft'52MASS? but the x- and y-axes were interchanged to 
avoid an infinite value of the slope (see Fig. 6). The resultant slopes Ak^/ Ex^^j^^^^Ks are 
also summarized in Table 1. 

The uncertainty of the RGB slopes is an error source of this method. Hence the fifth 
step described in §3 was carried out with different RGB slopes, which are la larger and 
smaller than those previously used. Here a was obtained when the histograms of (slope) ~^ 
are fitted with the Gaussian function. The changes of Aks/Eks-x are only less than a few 
% in the IRAC bands, although those in the 2MASS bands are 2 - 5%. 

The ratios of total to selective extinction Axg/ Ek^-x = Aks/^Aks — Ax) provide us 
with the ratios of absolute extinction Ak^ '■ Ax for the IRAC bands. Table 1 presents 
the extinction ratios Ax/Aks-i which are directly provided by Aks/Eks-x- We obtain the 
wavelength dependence of extinction between Ks and IRAC bands, Ajcg '■ ^[3.6] ■ ^[4.5] ■ 
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Fig. 5. — (a)-(d): Distribution of -ft's'smius magnitudes of a RC peak and -ft'ssiRius ~ ^ colors 
of RGB in the 0?2 x 0?2 fields for the [3.6] (a), [4.5] (b), [5.8] (c), and [8.0] (d) IRAC bands. 
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Fig. 6. — (a),(b): Distribution of -fi'ssiRius magnitudes of a RC peak and X — Kssimcs colors 
of RGB in the 0?2 x 0?2 fields for the 2MASS J (a), and H (b) bands, (c): Distribution of 
-^52MAss ~ -^5siRius colors of RGB and -ft'ssiRius magnitudes of a RC peak in the 0?2 x 0?2 
fields. To avoid an infinite value of a linear fitting, x- and y-axes were interchanged compared 
to (a) and (b). 
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A[5.8] ■■ ^[8.0] = 1 : 0.50 : 0.39 : 0.36 : 0.43. 

We also obtain the extinction ratio in the J, H and Kg bands in the 2MASS system, 
^J2MASS : ^i?2MASS ■ ^i^s2MAss = 2. 86 ± 0.08 : 1.60 ± 0.04 : 1. These are shghtly smaller 
than those obtained for the MKO system, Aj : Ah : Aks = 3.02 ± 0.04 : 1.73 ± 0.03 : 1 
(Nishiyama et al. 2006), but the differences of them are within less than 2a and 3a for 
^J2MASs/^Ks2MASs ^if2MASs/^Ks2MASS' respectively. 

To examine possible variations of the extinction law, we divided the survey area into 
quadrants, N+ (+3° > Z > 0°,+l° > 6 > 0°), S+ (+3° > Z > 0°,0° > 6 > -1°), N- 
(0° > Z > -3°, +1° > 6 > 0°), and S- {0° > I > -3°, 0° > 6 > -1°). The ratios Aks/Eks-x 
in the quadrants are listed in Tabic 2. Wc do not find significant deviation of the ratios 
from that for all the data points, but we find some trends. The ratios in N— tend to have 
a smaller value, while those in S— have a larger one. These trends, smaller N— and larger 
S— , have been obtained in Axg/ Eh-Ks (Nishiyama et al. 2006), and similar variations in 
the extinction law seem to be present among the quadrants in the MIR wavebands. We note 
again that we do not find statistically significant evidence for differing the extinction laws 
in different lines of sight toward the GC in the IRAC wavebands. 

Since most of stars we detected are in the bar structure whose major axis is oriented at 
~ 20° — 40° with respect to the Sun-Galactic center fine, different average distance of the 
giants in a given patch of the sky may cause systematic shifts in positions of the RC peaks 
and the RGB colors on the Ks versus Ks — X plot. As shown in Nishiyama et al. (2005), mean 
magnitude of RC stars weakly depends on the Galactic longitude. However, for changing 
the ratios of total to selective extinction Ak^/ Exg-x, it is required that the distance to the 
stars needs to be highly correlated with the reddening. Such correlation seems to be unlikely 
to exist^. A good test for the existence of this systematic error is to compare the ratios 
Axs/ Exg-x for the quadrants, because the error should be reduced in smaller regions due 
to smaller difference of the distance. As described above, we found insignificant deviation 
of Aks/Eks-x for the quadrants, and very small difference between the ratios for all the 
survey area and the weighted mean for the quadrants (see Table 2), suggesting very small 
systematic shifts on the Kg versus Ks — A plots due to the difference of the distance. 

To confirm that the current method using the RGB and RC stars in deriving the total 
to selective extinction ratios is consistent with the Nishiyama et al. (2006) method using 
only the RC stars, we re-analyzed the Nishiyama et al. (2006) data using the RGB and RC 



^We did not find any correlation between extinction values of the fields and their longitude at | ^ |< 2°. 
A weak correlation between them at | Z |> 2° can be found, but we confirmed that Aks/Eks-X for all the 
data points and that for | Z | < 2° are almost the same. 
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stars. The resultant ratios are Aks/Eh-Ks = 1-46 ± 0.03 and Aks/Ej^Ks = 0.499 ± 0.018, 
quite consistent with the previously derived ratios 1.44 ± 0.01, 0.494 ± 0.006 in Nishiyama et 
al. (2006). Therefore, the use of RGB colors does not affect the results. 

Table 1: The wavelength dependence of the interstellar extinction in the IRSF/SIRIUS 
{J,H,Ks\ similar to the MKO filter system), SST/lRkC ([3.6], [4.5], [5.8], [8.0]) and 2MASS 

(-^2MASS, -f^2MASS, -ft^S'2MASs) bauds. 



band 


wavelength [/xm] 


Aks/ Eks-x' 




J 


1.25 




3.02 ±0.04 


H 


1.63 ^ 




1.73 ±0.03 ^ 


Ks 


2.14 ^ 




1.0 


[3.6] 


3.545 ^ 


2.01 ±0.04 


0.50 ±0.01 


[4.5] 


4.442 d 


1.64 ±0.02 


0.39 ±0.01 


[5.8] 


5.675 ^ 


1.56 ±0.03 


0.36 ±0.01 


[8.0] 


7.760 


1.74 ± 0.04 


0.43 ±0.01 


•'^2MASS 


1.240 ^ 


-0.528 ±0.015 


2.89 ±0.08 


-f^2MASS 


1.664 


-1.61 ±0.04 


1.62 ±0.04 


-f^52MASS 


2.164 


-(0.012 ±0.002)-! 


1.012 ±0.002 



References. — ^ The ratio of total to selective extinction derived from the slopes in Fig. 5 and 6. ^ The 
relative extinction to the Ks band (IRSF/SIRIUS system). The mean effective wavelengths of the SIRIUS 
filters (see §5.1; Nishiyama et al. 2006). ^ The isophotal wavelengths of the 2MASS and IRAC filters (see 
Table 1; Indebetouw et al. 2005). ^ Prom Nishiyama et al. (2006). 



5. Discussion 

5.1. Comparison of Wavelength Dependence of Extinction with Previous 

Studies toweird the GC 

Fig. 7 shows the derived Ax/Ak^ toward the GC (this study; Rieke & Lebofsky 1985; 
Lutz 1999; Jiang et al. 2006). A simple power law of Ax oc A-^'^^ is also represented. Before 
the observations by using ISO/SWS, the wavelength dependence of interstellar extinction 
in the near- to mid- infrared wavelength range (1-6 /xm) was thought to be represented by 
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Table 2: The Axg/ Exg-x {Axg/ Ex^Ks ioi the 2MASS bands) for all the survey area, four 
quadrants, and their weighted mean 



band 


all data points^ 




S+b 




S-b 


weighted mean^ 


[3.6] 


2.01 ±0.04 


2.01 ±0.07 


2.03 ±0.09 


1.85 ±0.06 


2.05 ±0.09 


1.96 ±0.04 


[4.5] 


1.64 ±0.02 


1.62 ±0.04 


1.60 ±0.05 


1.59 ±0.06 


1.69 ±0.07 


1.62 ±0.02 


[5.8] 


1.56 ±0.03 


1.56 ±0.04 


1.52 ±0.05 


1.44 ±0.04 


1.64 ±0.07 


1.52 ±0.03 


[8.0] 


1.74 ± 0.04 


1.82 ±0.06 


1.67 ±0.07 


1.80 ±0.08 


1.83 ±0.09 


1.77 ±0.04 


J2M 


0.528 ±0.015 


0.497 ±0.029 


0.534 ± 0.065 


0.494 ± 0.021 


0.587 ±0.041 


0.510 ±0.016 


H2M 


1.61 ± 0.04 


1.64 ±0.05 


1.54 ±0.08 


1.48 ±0.07 


1.63 ±0.10 


1.58 ±0.04 



References. — ^ The same value listed in Table 1. ^ Aks/Eks-\ for ~ 3°x ~ 1?5 quadrants (see §4). 
Weighted means of the values for N+, S+, N— , and S— . 

a power law and to be "universal" (Draine 1989). However, the observation of HII regions 
around Sgr A* by using ISO/SWS shed serious doubt on a simple power law, and recent 
observations toward star forming regions with SST/IKAC show clear discontinuity of the 
power law, preferring flatter extinction toward longer wavelengths. 

Analysis of the spectrum of the GC obtained with ISO/SWS revealed an extinction law 
characterized by a relatively flat behavior at 3 - 8 //m (Lutz et al. 1996). The extinction 
measurements were improved by Lutz (1999, open diamonds in Fig. 7), which reinforce the 
previous flat extinction. In Fig. 7, their A\/Av is converted to Ax/Ak^ by assuming a 
_^-i.99 extinction law (Nishiyama et al. 2006) from the Ks wavelength to their A = 2.625 /xm 
point. Fig. 7 shows that the extinction law derived by Lutz (1999) is very similar to those 
presented in this work, but discrepancy is slightly large in the [4.5] and [5.8] bands. Some 
of the discrepancies in ~ 3 — 6 //m wavelength between Lutz (1999) and our results may 
be explained by large absorption features observed toward the Sgr A*. Chiar et al. (2000) 
analysed spectra (2.4 - 13 jim) toward the Sgr A* and two sources in the Quintuplet cluster, 
and discussed composition of dust along the hues of sight to them. We can see a deep 
and broad absorption feature around ~ 6 ^m. toward the GC for which Chiar ct al. (2000) 
proposed ice mixtures (H20:NH3:C02) and HCOOH. A deep absorption feature of H2O ice 
at 3 jim. was also detected. The same spectrum was used for Lutz (1999), and thus their 
higher extinction at 3.0392, 3.2970, and 5.9082 fim. can be attributed to the ice along the line 
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of sight. In contrast, a relatively shallow and narrow absorption was observed for the sources 
in the Quintuplet cluster. Although it is still unclear how deep and broad the features are 
along the other lines of sight in the GC region, and two of the higher data points (4.3765 
and 5.1287 /im) cannot be attributed to such absorption features, if the deep and broad 
absorption features at ~ 3 and 6 /zm are characteristic to only Sgr A* and its immediate 
vicinity, some of the discrepancies can be naturally explained. 

Fig. 7. — The wavelength 
dependence of extinction 
(Ax/Aks) toward the Galac- 
tic center region. The ratios 
derived by this study {filled 
blue circles), Lutz (1999, open 
diamonds), Jiang et al. (2006, 
open triangle), and Rieke & 
Lebofsky (1985, filled squares) 
are represented. Note that in 
the case of Rieke & Lebofsky 
(1985), Ax/Ak instead of Ks 
is represented. A simple power 
law. Ax oc (Draine 1989) 

is shown by the dashed line. 
Although it is not toward the 
Galactic center, the extinction 
curve derived by Indebetouw et 
al. (2005, equation 4) is shown 
by thin red smooth line for 
reference. 

Fig. 7 shows that the extinction law derived by Rieke & Lebofsky (1985) is much smaller 
than Lutz (1999) and our results at > 4 /xm, and the decrease in the 1 - 2.5 /im range is also 
very different from our results. Cardelli et al. (1989) fit the data of Rieke & Lebofsky (1985) 
with a power law A~^'^^, and it has quite often been referred to as the standard extinction 
law. However, since Rieke & Lebofsky (1985) laid great emphasis on the determination of 
Av/Eb-v in accordance with the reddening data taken before, J magnitudes of o Sco, Cyg 
0B2 No. 12, and only two stars near the Galactic center were measured with relatively large 
uncertainties in the color excesses mainly due to the assumed errors in intrinsic colors of 
these stars. Therefore, the color excesses they determined Ey-j / Eb~v = 2.19 ± 0.04 and 
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Ev-h/Eb-v — 2.55 ± 0.02 on the assumption of Ev-k/Eb-v — 2.744 might be in fact 
compatible with the current results if we allow for these uncertainties. 

The extinction ratios toward the GC in the J, H and Ks bands in the 2MASS system 
are consistent with the IRSF/SIRIUS results, as shown in Table 1. These ratios are plotted 
in Fig. 8, and are consistent with a power law Ax oc A~^ °, together with the Nishiyama 
et al. (2006) results. Therefore, it is established that the interstellar extinction toward the 
GC in the range of J, H and Ks bands is well fit by a steep power law of A"^ *^, not A^^'^. 
Of course, the power law is only an approximation. Since 1) the extinction shows positive 
deviation from the power law as the wavelength increases to > 3 fim. as we sec above, and 2) 
a number of dust grain models (e.g., Weingartner & Draine 2001; Dwek 2004; Voshchinnikov 
et al. 2006) predict such positive deviation, it is possible that the power index changes in 
the J, H, and Ks wavelength range. 



3.0 



^ 2.0^ 



1.0 



•Nishiyama et al. 2006 

(IRSF/SIRIUS) 
OThis work (2MASS) 




Fig. 8. — The wavelength 
dependence of extinction 
(Ax/Aks) in the J, H, and 
Ks bands toward the GG. 
The results of Nishiyama et al. 
(2006) in the MKO system are 
presented by filled circles, and 
that in the 2MASS system by 
open circles. A simple power 
law Ax (X A-2°o fitted to all 
these points is shown by the 
straight line. Horizontal bars 
indicate the filter widths. 



Wavelength [fim] 



5.2. Compcirison in the Ratios of Color Excesses 



A number of authors obtained the ratios of color excesses first, and then derived the 
wavelength dependence Ax/Axg of the interstellar extinction. However, the assumed ratio 
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-^hI-^Ks leads to a quite different wavelengtfi dependence, as is evident from the equation 
(2) and was already pointed out by Nishiyama et al. (2006) and Flaherty et al. (2007). 
Therefore, we convert our A\/ A^g to the ratios of color excesses and compare them with the 
previous results (Fig. 9). This process causes larger errors because of the error propagation, 
but not any uncertainty at all due to assumption of an unknown parameter as in the inverse 
procedure. 

Flaherty ct al. (2007) and Roman-Ziiriiga et al. (2007) obtained E\_Ks/ ^h-Ks (o^ its 
reciprocal) from the slopes of the distributions of stars in color-color diagrams. Indebetouw 
et al. (2005) measured the color excess ratios E\^Ks/ Ej-k^^ E\-Ks/ Eh-Ks was derived 
from color-color diagrams [open squares) by Flaherty et al. (2007, see their Table 2) using 
the same data sets of Indebetouw et al. (2005) and the same method as used in Flaherty et 
al. (2007). Thus, direct comparison of the wavelength dependence of extinction is possible 
in color-excess ratios. In Fig. 9, the color-excess ratios Ex-Ks/Eh-Ks are plotted against 

following the convention to compare with scattering calculation. 

Differences in E^-Ks/ Eh-Ks seem to exist between the star forming regions and the 
diffuse interstellar medium. The ratios toward the star forming dense core (Roman-Ziiniga 
et al. 2007) and the nearby star forming regions (Flaherty et al. 2007) are very similar, and 
are in very good agreement within their errors. As indicated by Flaherty ct al. (2007), we 
can find a systematic separation of the ratio toward the off-cloud regions (Indebetouw ct al. 
2005) from toward the star forming regions. The ratios obtained in this study (calculated 
from Ax/Aks) tend to have lower values than those for star forming regions, in the three 
bands but [8.0]. The higher extinction in the [8.0] band in the direction of the GC than 
that toward I — 284° (Indebetouw et al. 2005) is consistent with the claim that the GC line 
of sight has larger silicate absorption relative to Ay (larger Tg.r/Ay; Roche & Aitken 1985) 
than other hues of sight (see also Table 1, Draine 2003). A separation between the extinction 
laws for molecular clouds and for the diffuse interstellar medium in the color excess ratios 
-E'[3.6]-[4.5]/-E'[4.5]-[5.8] Suggested by Flaherty et al. (2007) is not confirmed in Fig. 9 due to the 
large uncertainties. However, we can find a clear difference in Ex-KsI Er-Ks between the 
star forming regions and the off-cloud regions including toward the GC. 



5.3. Absolute Extinction Ratios in Different Lines of Sight 

Absolute extinction has only been measured for the GC, and comparison with other lines 
of sight requires an assumption of at least one extinction ratio (e.g., Ah/Ak^)- Here, we 
assume a single power law for the interstellar extinction instead, and try to derive its index. 
Flaherty et al. (2007) measured Ej_Ks/Eh-Ks in NGC 2024/2023, NGC 2068/2071, and 
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1 / 1 [|im]"^ 

Fig. 9. — Color excess ratio Ex^Ks/ ^h-Ks versus wavelength diagram in the IRAC bands 
derived by this study {blue filled circles) and Roman- Ziifiiga et al. (2007, open circles). The 
ratios derived by Flaherty et al. (2007) for five star-forming regions {solid lines) and their 
weighted mean {filled triangles) are also plotted. Here all the H and Kg bands are those of 
the 2MASS survey. The ratio Ex-Ks/ ^h-Ks '^'^^ directly measured by Indebetouw et al. 
(2005), but Flaherty et al. (2007) directly determined Ex-KsI Er-Ks for the same data sets 
of Indebetouw et al. (2005), and this is shown by red open squares. So the data points shown 
in this figure, except ours (blue circles), are directly measured from color-color diagrams. 
Only ours are converted from Ax/Axg to Ex-Ks/^h-Ks- 



-22- 



Serpens, which are 3.00 ± 0.04, 2.98 ± 0.03, and 3.05 ± 0.04, respectively. They also derived 
the ratio for the / = 284° region (Indebetouw et al. 2005) with the result of 3.07±0.04, point- 
ing that no significant variation seems to exist in the near-infrared extinction law through 
these molecular clouds and through the diffuse interstellar medium. Assuming the effective 
wavelengths of the 2MASS observations to be f.240, f .664, and 2.f64 fim (Indebetouw et 
al. 2005) and that the extinction is approximated by a power law, we can derive the power 
law index for these lines of sight. The resultant indices are from —2.04 to —2.18. Thus, the 
interstellar extinction for all these hues of sight seems to be fitted by a single power law of 
Ax oc A~^ °, which is the same as that in the GC, very well in the J, H, and Ks wavelength 
range. However, the extinction law derived by Indebetouw et al. (2005) decreases much more 
slowly to the longer wavelength (slower than A~^'^), and nonetheless their data are shown 
by Flaherty et al. (2007) to have a consistent Ej^KsI ^h-Ks ratio 3.07 ± 0.04. (Note that 
Table 1 in Indebetouw et al. (2005) lists smaller Ej^KsI ^h-Ks ratios, possibly because of 
the difference in the selection criteria of background stars by them and by Flaherty et al. 
(2007).) Therefore, from the color excess ratios one cannot resolve the degeneracy. Both a 
steep power law upto the Kg band similar to the GC and a much gentler extinction curve 
are consistent with the color excess ratios. In particular, the frequently employed JHKs 
color-color diagrams can be quite insensitive to the change in extinction law. Hence, we 
cannot determine a value of Ah/Aks even from a seemingly similar Ej^KsI Eh-Ks values. 
We further should be cautious, because line-of-sight differences in Ej_Kg/ Eh-Ks do seem 
to exist. For example, Naoi et al. (2006) found color excess ratios in p Oph and Cha I 
star-forming regions consistent with that in the GC, but a significantly larger ratio in the 
Coalsack Globule 2 (Naoi et al. 2007). 

It is quite straightforward to determine A\/ Ak^ in this work from the direct measure- 
ment of the ratios of total to selective extinction A^sl Eks-\- In contrast, to derive Ax/A^s 
in the equation (2) in the color-color method is a difficult and crucial task, which Indebe- 
touw et al. (2005) made from the determination of An/ Axg. Since Flaherty et al. (2007) and 
Roman-Ziiiiiga et al. (2007) did not determine An/AKg-, they converted their color excess 
ratios to Ax/Ak^ using the ratio Ah/Aks — 1-55 derived by Indebetouw et al. (2005). How- 
ever, there is no clear evidence for universality of Ah/Akqi and a different Ah/Aks leads to 
a different Ax/Axg as described in §5.2, and as also pointed out by Flaherty et al. (2007). 
The difference is shown in Fig. 10 in which Ah/Aks = 1.60 in addition to 1.55 is used to 
derive another set of Ax/Ak^ ratios based on the Flaherty et al. (2007) and Roman-Ziiniga 
et al. (2007) color excess observations. 

The ratio of Ah/Ak^ = 1.55 produces larger separation between the line of sights to 
molecular cloud and diffuse interstellar medium. The average of the four lines of sight 
to the star-forming regions studied by Flaherty et al. (2007) and the dark cloud core by 
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Roman-Zuniga et al. (2007) show higher extinction in the IRAC bands than the GC and 
the / = 284° diffuse medium (Indebetouw et al. 2005). In contrast, if we assume the ratio of 
An/Axg = 1.60, the molecular-cloud data points shift to the GC and the diffuse interstellar 
medium, although An/Axg ~ 1.69 is needed to make the results of Flaherty et al. (2007) 
and Roman-Ziihiga et al. (2007) agree with our results. Therefore, if we assume that the 
near-infrared extinction laws are different in the dense clouds and diffuse clouds, then the 
difference continues to the mid-infrared, and vice versa. Currently, we do not have sufficient 
data to distinguish these two possibilities because the measurement of the total to selective 
extinction is very difficult in the hues of sight other than the GC. However, the difference 
in Ex-Ksl Er-Ks (see Fig. 9) and the required Ah/Aks of ^ 1.69 suggest that all of the 
observational results discussed in this paper cannot be reconciled with a single extinction 
law. 



6. Conclusion 

We have measured the wavelength dependence of interstellar extinction toward the GC 
in the 1.2 - 8.0 /xm region by combining the IRSF/SIRIUS infrared observations and the 
2MASS and SST/IKKC catalogs. The extinction in the wavelength range of J, if, and Kg 
is well fitted by a power law of steep decrease A\ oc A~^''^ toward the GC. Furthermore, 
the fiattcning of the extinction from a simple extrapolation toward the longer wavelength of 
the power law at < 3 iim has been confirmed. In particular, the extinction has a relatively 
shallow and broad minimum in the SST/IUKC wavelength range; 5] and ^[5.8] are only 
slightly smaller than 0.4 times Ak^- 

This dependence has been derived directly from the observation of reddening in pro- 
portion to the Kg extinction for the first time involving the SST/IR.AC wavebands. In this, 
the Ks magnitudes of RC stars and the Ks — A colors of RGB stars serve as a tracer of 
the reddening vector in the color-magnitude diagrams, and the ratios of total to selective 
extinction Axg/ Eks-\ have been obtained by a variant of the "RC method" originated from 
the variable-extinction method. It is interesting to note that this method is only sensitive 
to the spatially variable component in the surveyed area; when a star cluster suffers patchy 
0.8 — 1.4mag extinction from region to region, then the total to selective extinction ratio for 
the variable 0.6mag component is derived, but the characteristics of the ubiquitous 0.8 mag 
extinction component remain unknown. This can be its strong point because it is not so 
sensitive to the intrinsic magnitudes and colors of these background stars; on the other hand, 
the location of the dust grains causing this variable extinction in the long line of sight to the 
GC is not certain. However, since the silicate absorption seems fairly strong, the dominant 
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Fig. 10. — Left : The wavelength dependence of extinction (Ax/Aks)- The blue filled circles 
(this study) is the extinction law toward the Galactic center, and the red open squares are that 
for off-cloud regions in the Galactic plane (Indebetouw et al. 2005). Those for star forming 
regions derived by Flaherty et al. (2007) and Roman-Zuniga et al. (2007) are shown by filled 
triangles and open circles, respectively, in the assumption of Ah/Aks = 1-55 (Indebetouw et 
al. 2005). A simple power law. Ax (x A~^'^^ is shown by the dashed line (Draine 1989). Right 
: The same extinction laws derived by this work and Indebetouw et al. (2005) are shown, 
but those for Flaherty et al. (2007) and Roman-Ziiniga et al. (2007) are slightly different due 
to the different assumption of Ah/Ak^ = 1.60 corresponding to a steep power law. 
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grain population probed in this method is probably in the vicinity of the GC. It should be 

also noted that the assumptions used here are that the RC stars and RGB stars have spatial 
distributions with their centers in common, and that the reddening of RGB stars can be 
measured precisely from the shape of upper giant branch. 

The wavelength dependence seems to be different among the various lines of sight. 
In particular, the GC and off-cloud regions show different extinction law from star-forming 
regions, although the difference is not very large compared to the observational uncertainties. 
The general behavior of extinction curves, which decreases steeply in the near-infrared, but 
decreases only slightly beyond the Ks band, is an important constraint when dust grain 
models are discussed. 
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